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We evaluate the shadowing effect in deeply virtual and real Compton scattering on nuclei in the 
framework of the color dipole model. We rely on the soft photon wave function derived in the 
instanton vacuum model, and employ the impact parameter dependent phenomenological elastic 
dipole amplitude. Both the effects of quark and the gluon shadowing are taken into account. 
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Compton scattering, 7* + p — > 7 + p, with initial photons both real or virtual, has been a subject of intensive 
theoretical and experimental investigation [l|-[l6j. While in the case of deeply- virtual Compton scattering (DVCS), 
where the initial photon is highly- virtual, the QCD factorization has been proven [j| 0, H| and the amplitude can be 
expressed in terms of the generalized parton distributions (GPD) [H-[l5j|. in the case of real Compton scattering (RCS) 
■ the available theoretical tools are rather undeveloped. 

On the one hand, as it has been shown in jl7l lisj j , for large momentum transfer A_l it is possible to factorize the 
RCS amplitude [19L |2G | and express it in terms of the distribution amplitudes of the proton. On the other hand, it is 
possible to express the amplitude of the process via the minus 1st- moment of GPDs at zero skewedness [f| [2l|, [22[ ■ 

DVCS and RCS on a proton have been studied recently within the color dipole approach in [2^, [24j ■ Here we extend 
that study to nuclear targets. The DVCS process on a nuclear target has been measured at HERA by HERMES 
collaboration (2(| and may be also studied at the future Electron Ion Collider (EIC) and Large Hadron Electron 
Collider (LHeC) [271 l28l|. The Real Compton Scattering (RCS) may be measured at the LHC, as a subprocess in 
hadron-hadron collisions in ultraperipheral kinematics. Since both DVCS and RCS are studied in the high-energy 
kinematics, the nuclear effects reveal themselves as shadowing corrections. 

The general framework for evaluation of the shadowing corrections is the Gribov-Glauber approach [29(. While in 
the asymptotic high-energy ( "frozen" ) regime the shadowing corrections were studied in [2^, [30] , we use an approach 
which is also valid for intermediate energies. Also, we take into account the gluon shadowing corrections, which appear 
f° r xb ^ 1CP 3 and give a sizeable contribution for xb ~ 10~ 5 . 

The paper is organized as follows. In Sections [TTI we review the general formalism of the color dipole approach. In 
, Section Hill we discuss the frequently used frozen approximation which is valid for asymptotically large energies. In 
Section [TV] we discuss the method which will be used for calculations of nuclear shadowing effects and demonstrate 
that for asymptotically large energies it reproduces the results from Section IIIII In Section [V] we discuss the gluon 
shadowing and its effect on the DVCS and RCS observables. In Section I VII the wavefunction of a real photon is 
evaluated in the instanton vacuum model. In Section I VIII we present the results of numerical evaluation, and in 
Section [Villi we draw conclusions. 



II. COLOR DIPOLE MODEL 



The color dipole model is particularly efficient at high energies, where the dominant contribution to the Compton 
amplitude comes from gluonic exchanges. Then the general expression for the Compton amplitude on a nucleon has 
the form, 



A$> (s, A) = e«e« / dpxdfotfricPr^f (ft, r 2 ) A d (ft, ff; fa r 2 ; A) ^ (ft, n) , 



(1) 
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where e$ is the photon polarization vector; ft : 2 are the light-cone fractional momenta of the quark and antiquark, f'l 2 
are the transverse distances in the final and initial dipoles respectively; A is the momentum transfer in the Compton 
scattering, A d (...) is the scattering amplitude of the dipole on the target (proton or nucleus), and ^^(n (ft r) are the 
wavefunctions of the initial and final photons in the polarization state i (3lj . 

At high energies in the small angle approximation, A/y/s <C 1, the quark separation and fractional momenta j3 are 
preserved, so, 
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is the gluon GPD of the target, P' = P + A,_P = (P + P')/2, G^ v (x) is the gluon loop operator, £qq (x, y) is the 
Wilson factor required by gauge covariance. For this GPD we use a gaussian parameterization [55l l56l |66|. 




where o'o( :c )j -RoW' are the free parameters fixed from the DIS and irp scattering data. We shall discuss them 

in more detail in Section [VIII The parameterization (J5J) does not depend on the longitudinal momentum transfer and 
decreases exponentially as a function of A 2 . Since the parameterization ([5]) is an effective one and is valid only in the 
small- x region, we do not assume that it satisfies general requirements, such as positivity [65| and polynomiality Q 
constraints. 

The prefactor (fc + —) (fc — —) in ([3]) guarantees convergence of the integrals in the parameterization @. In 
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the forward limit, the amplitude reduces to the saturated parameterization of the dipole amplitude proposed by 
Golec-Biernat and Wusthoff (GBW) @, 

*,(/9,r) = 2y'd 2 6'3m/,>^/3) = ^y'^ C ^ (6) 

Generally, the amplitude /^(■••) involves nonperturbative physics, but the asymptotic behaviour for small r is con- 
trolled by pQCD (32: 

fg{r,A,p) r ^<xr\ 

up to slowly varying corrections ~ ln(r). 

Calculation of the differential cross section also involves the real part of scattering amplitude, whose relation to the 
imaginary part is quite straightforward. According to [33], if the limit lim ( Im J ) exists and is finite, then the real 
and imaginary parts of the forward amplitude are related as 



Kef (A = 0) = s Q tan 
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In the model under consideration the imaginary part of the forward dipole amplitude indeed has a power dependence 
on energy, Imf(A = 0; s) ~ s Q , so ([7]) simplifies to 

^=fn( f ( Q -l)) S «. (8) 
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This fixes the phase of the forward Compton amplitude, which we retain for nonzero momentum transfers, assuming 
similar dependences for the real and imaginary parts. Finally we arrive at, 



^ = {e + i)ef{q')e^{q) I d 2 r / dfivftf, r)*V>(J3, r) 3m/^(f, A, ft s), 



,00, 



(9) 



For the cross-section of unpolarized Compton scattering, from (j9]) we obtain, 
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III. NUCLEAR SHADOWING IN THE "FROZEN" LIMIT 



Nuclear shadowing signals the closeness of the unitarity limit. Hard reactions possess this feature only if they have 
a contribution from soft interactions. In DIS and DVCS the soft contribution arises from the so called aligned jet 
configurations [34], corresponding to qq fluctuations very asymmetric in sharing the photon momentum, /3 <C 1- Such 
virtual photon fluctuations, having large transverse separation, are the source of shadowing [40(1 . 

Calculation of nuclear shadowing simplifies considerably in the case of long coherence length |35| . i.e. long lifetime of 
the photon fluctuations, when it considerably exceeds the nuclear size. In this case Lorentz time dilation "freezes" the 
transverse size of the fluctuation during propagation though the nucleus. Then the Compton amplitude of coherent 
scattering, which leaves the nucleus intact, has the same form as Eq. © with a replacement of the nucleon Compton 
amplitude by the nuclear one, 
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where b is impact parameter of the photon- nucleus collision, T A (b) = dz p A (b, z) is the nuclear thickness function, 
given by the integral of nuclear density along the direction of the collisions. In this expression we neglect the real 
part of the amplitude which is particularly small for a coherent nuclear interaction. 

For incoherent Compton scattering, which results in nuclear fragmentation without particle production (quasielastic 
scattering), the cross section has the form [3r| . 
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Here B e i is the i-slope of elastic dipole-nucleon amplitude. In this equation we treated the term quadratic in the 
dipole cross section as a small number and expanded the exponential in curly brackets. 



IV. ONSET OF NUCLEAR SHADOWING 
A. Coherent Compton scattering 



The regime of frozen dipole size discussed in the previous section is valid only at very small xb in DVCS, or at high 
energies in RCS. However, at medium small xb a dipole can "breath", i.e. vary its size, during propagation through 
the nucleus, and one should rely on a more sophisticated approach. 
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In this paper we employ the description of the onset of shadowing developed in [37[ and based on the light-cone 
Green function technique (38[. The propagation of a color dipole in a nuclear medium is described as motion in an 
absorptive potential, i.e. 



dW (z 2 ,r 2 ;zi,ri) A r2 W (z 2 ,r 2 ; z 1 ,r 1 ) _ ip A (z 2 , r 2 ) <7g q (r 2 ) 
dz 2 va(l — a) 2 



W{z 2 , r 2 ;zi,r 1 ), (13) 



where the Green function W (z 2 , r 2 \ z\, r\) describes the probability amplitude for propagation of dipole state with 
size r\ at the light-cone starting point z\ to the dipole state with size r 2 at the light-cone point z 2 . Then the shadowing 
correction to the amplitude has the form 
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Equation (|13[) is quite complicated and in the general case may be solved only numerically |39j |. However in some 
cases an analytic solution is possible. For example, in the limit of long coherence length, l c 3> Ra, relevant for 
high-energy accelerators like LHC, one can neglect the "kinetic" term oc A r2 W (z 2 , r 2 ; Z\, rj.) in (fT5]l . and get the 
Green function in the "frozen" approximation [38j, 



W (z 2 ,r 2 ; Zl ,n) = S 2 (r 2 - n) exp ^-^<Tq q (n)J d(p A ((,b) | . (14) 
Then the shadowing correction ([T5|) simplifies to 
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(a,r)exp (-^ 9q (r) Jd(p A ((, b) j % n (a, r) e *mi«(«-*i>. 

If we neglect the real part of the amplitude and the longitudinal momentum transfer k m in (which is justified for 
asymptotically large s), and average over polarizations, then taking the integral over zi >2 "by parts" in (I15p . we get 
for the elastic amplitude 
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A(s,A x ) = 2 j d 2 be lK ^ J dad 2 r^ f (a,r) 1 - exp (r) J d(p A ((,b) 

L \ -oo 
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Another case where an analytical solution is possible is when the effective dipole sizes are small and the function 
Uq q (r) may be approximated as 

aq q {r) « Cr 2 . (17) 

This approximation cannot be precise even at high virtuality Q 2 in DVCS, since there are contributions of the 
aligned jet configurations mentioned above, which permit large dipoles even for lar ge Q 2 . Moreover, such aligned 
jet configurations of the dipole provide the main contribution to nuclear shadowing 40J. Nevertheless, for the sake 
of simplicity we rely on this approximation in order to estimate the magnitude of the shadowing corrections. Then 
Eq. (|13[) yields for W (z 2 , r 2 ; z%, r\) the well-known evolution operator of harmonic oscillator, although with complex 
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frequency 
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is comparable with the mean inter-nucleon spacing in nuclei and is much smaller than the radii of heavy nuclei. 
Therefore the frozen approximation employed in [25| cannot be used for an interpretation of HERMES data, instead 
one should rely on the Green function method described above. 



B. Incoherent scattering 

In addition to the coherent processes which leave the recoil nucleus intact, a large contribution to photon production 
comes from incoherent Compton scattering, where the target nucleus breaks up. Using the missing mass technique 
one can select events jA — > "/A* in which the nucleus breaks up to fragments without production of mesons. In this 
case one can employ completeness of the final states, which greatly simplifies the calculations. 

The analysis of such processes in electroproduction of vector mesons was done in 41[, and its extension to the 
DVCS and RCS is quite straightforward and yields 
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At sufficiently high energies one can rely on the frozen approximation introduced in the previous section, and this 
formula may be simplified [42| . 
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Correspondingly the cross section (fl9|) takes the form, 
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This expression reproduces Eq. (fl2j) derived in the "frozen" limit. It is easy to see that in the limit of a transparent 
nucleus, <Jq q T A "C 1, the cross section Eq. (j2"4l rises linearly with A. However in the limit of a very opaque nucleus 
(black disk limit), a qq T A ^> 1, the absorptive exponential factor in (|24|) terminates the contribution of central impact 
parameters, and da / dt\A ± =o oc A 1 / 3 . 

In case of the approximation (|17[) . we may use the explicit expression (|18[) for the Green function W (z 2 , f 2 \ zi, f\) 
to get 
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GLUON SHADOWING 



It has been known since |43| that in addition to the quark shadowing inside nuclei there is also a shadowing of 
gluons, which leads to attenuation of the gluon parton distributions. While nuclear shadowing of quarks is directly 
measured in DIS, the shadowing of gluons is poorly known from data @|45j], mainly due to the relatively large error 
bars in the nuclear structure functions and their weak dependence on the gluon distributions, which only comes via 
evolution. The theoretical predictions for the gluon shadowing strongly depend on the implemented model-while for 
x b ^ 10~ 3 they all predict that the gluon shadowing is small, for xb < 10 -3 the predictions vary in a wide range 
(see the review [45[ and references therein). Since in this paper we also make predictions for the LHC energy range, 
the gluon shadowing corrections should be taken into account as well. 

In the framework of the color dipole model the gluon attenuation factor R g was evaluated in the Gribov-Glauber 
approach in |46| . It was found convenient to evaluate the gluon attenuation ratio R g defined as 



T A {b)G N (x,Q 2 ,6)' 



where Gjv (x, Q 2 ,&) is the impact parameter dependent gluon GPD, relating it to the shadowing corrections in DIS 
with longitudinally polarized photons, 
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where Acr^ P = < J \ A ~ Aa]^ p is the shadowing correction at impact parameter b, and al p (x, Q 2 ) is the total 
photoabsorption cross section for a longitudinal photon. The process with longitudinal photons is chosen because the 
aligned jets configurations are suppressed by powers of Q 2 , so that the average size of the dipole is small, (r 2 ) ~ l/Q 2 , 
and nuclear shadowin g m ainly originates from gluons. 
As it was shown in |46| , 
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where p A {b, z) is the nuclear density, and T (x, Q , Azj is defined as 
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For heavy nuclei we may rely on the hard sphere approximation, p^t(r) ~ Pa (0)6 (Ra ~ T \ an d simplify 
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where L = 2^J R 2 A — b 2 . For the total cross-section after integration over J d 2 b we may get 
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The results of evaluation of the gluon shadowing are presented in Section IVII 
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VI. WAVEFUNCTIONS FROM THE INSTANTON VACUUM 

In this section we present briefly some details of evaluation of the wavefunction in the instanton vacuum 
model (see [471449] ] and references therein). The central object of the model is the effective action for the light 
quarks in the instanton vacuum, which in the leading order in N c has the form [48l 149} 

S = J d 4 x Cy In A + 2$ 2 (x) — ip(p + v — m— cLf ® $ • T m ® fL) i/^ , 

where T m is one of the matrices, T m — 1, ir, 75, 1x75, ip and $ are the fields of constituent quarks and mesons 
respectively, N/V is the density of the instanton gas, v = is the external vector current corresponding to the 

photon, L is the gauge factor, 



L(x,z) = Pexp i / (0 



(29) 



which provides the gauge covariance of the action, and f(p) is the Fourier transform of the zero-mode profile. 

In the leading order in N c , we have the same Feynman rules as in perturbative theory, but with a momentum- 
dependent quark mass /i(p) in the quark propagator 
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The mass of the constituent quark has a form 

(i(p)=m + Mf 2 (p), 

where m s=s 5 MeV is the current quark mass, M ~ 350 MeV is the dynamical mass generated by the interaction with 
the instanton vacuum background. Due to the presence of the instantons the coupling of a vector current to a quark 
is also modified, 

v = -> V = v + V nonl , 

ynonl w _ 2Mm ^M v ^ q) + ( g 2y (31) 

Notice that for an arbitrary photon momentum q the expression for V nonl depends on the choice of the path in l|29p 
and as a result one can find in the literature different expressions used for evaluations [3l|, [50l - [5^ | . In the limit p — > oo 
the function f(p) falls off as ~ so for large p ^> p^ 1 , where p ps (600 MeV)^ 1 is the mean instanton size, the 

mass of the quark p,{p) w m and vector current interaction vertex V ~ v. However we would like to emphasize that 
the wavefunction ^f(/3, r) gets contribution from both the soft and the hard parts, so even in the large-Q limit the 
instanton vacuum function is different from the well-known perturbative result. 

We have to evaluate the wavefunctions associated with the following matrix elements: 
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where T is one of the matrices T = 7^, 7^75, er^. In the leading order in N c one can easily obtain 
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e^<5 (V -(P + l) 1 + ) Tr (S(P)VS( P + q)T) . (33) 



The evaluation of (I3"c>l) is quite tedious but straightforward. Details of this evaluation may be found in [3~i| . 
The overlap of the initial and final photon wavefunctions in (jlpp was evaluated according to 

r, g 2 = 0)M>« (/3, r, Q 2 ) = £ ij (/3, r*, 0) J r (/3, r, Q 2 ) , (34) 

r 

where summation is done over possible polarization states T — '■y^, 'y^'ys; cr^y. In the final state we should use 
r* = r M + n^^J-Ii = r M — n M A g + r ~ L ; which is related to the reference frame with = in which the components ([55)1 
were evaluated. 

VII. NUMERICAL RESULTS 

In this section we present the results of numerical calculations. In this paper we consider two processes-DVCS 
and RCS. While physically they differ only by the kinematics, the parameterizations used for scattering with soft 
and hard photons are different. For example, for DVCS we have photons with large Q 2 -in this region we have 
Bjorken scaling, so all the model parameters such as basic cross-section ctq and saturation radius Ro in Eq. 
should depend on the Bjorken xb- A widely accepted parameterization which incorporates this features is the GBW 
parameterization (53| . On the contrary, for RCS when Q 2 is vanishingly small, xb vanishes and thus all variables 
should depend on Mandelstam s. An example of such parameterization is KST parameterization introduced in [54T4571 ] . 
Up to the best of our knowledge there is no single parameterization which incorporates both asymptotic cases. For 
this reason for DVCS we used the GBW-style parameterization [5^, [66| for the nonintegrated gluon density Eq. (J5J) , 
which has a form 

(Jo{x) — 23.03 mb = const, (35) 
R (x) = OAfmx (x/x )°- 1A4 , (36) 

B(x) = B Tp ^ pp - i R 2 (x), (37) 
where x = 3.04 x 10~ 4 , and B^ p ^ pp (x,Q 2 > lGeV 2 ) « 5GW" 2 H. 
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Figure 1: The nucleus to nucleon cross section ratio for coherent DVCS as function of different kinematical variables. Left: 
^-dependence of the shadowing, t = t m in, for different Q 2 and A. From bottom to top: Q 2 — 5 GeV 2 and Q 2 = 20 GeV 2 . 
Right: Q 2 -dependence, t — — O.OlGeV 2 , for different Xb and A. From bottom to top: xb = 10 -5 and xb = 10~ 3 . 



For RCS we used a KST-style parameterization, which has the form [23j, |24|, I54H571 16 



23.6 



•so 



Ro(s) = 0.88 fm 



8 rl 

0.079 



so 



1.45 



(?) 



0.45 



s + si 
1 



(38) 

(39) 

(40) 
(41) 



where sq ~ 1000 GeV 2 , s± w 3600 GeV 2 . Note that in the large-s region considered in this paper we may neglect 



the second (Reggeon) term in Eq. (|39l) and set si « in (|40|) . As one can see from the Figure [1] the shadowing 
correction is increasing towards small xb, and for xb ^ 10~ 5 the nuclear cross-section ratio decreases by a factor of 
two compared to the naive estimate da a ~ F\(t)doN- As a function of the momentum transfer t, the shadowing 
correction reveals the behavior qualitatively similar to the nuclear formfactor FA(t): it steeply drops at small-i and 
has zeros for some t. Notice, however, that the zero positions in the cross-section do not coincide with the zeros of 
the formfactor. This is a result of shadowing which suppresses the contribution of the central part of the nucleus and 
modifies the 6-dependence of the cross section compared to the formfactor. Notice that, as was discussed in [59j . if we 
had sufficiently high resolution in t, it would be possible to measure the contributions of the "pure" DVCS, without 
contributions from the Bethe-Heitler, near the zeros of the formfactor. 

The Q 2 -dependence of the nuclear ratio is shown in the right pane of the Figure [T] 

In the Figure [5] we compare the results for the coherent cross section ratio evaluated with and without gluon 
shadowing. As expected, gluon shadowing is very small at xb > 10~ 3 and gives ~ 20 — 30% contribution for 
x ~ 10~ 5 . Similar dependence is observed for all the other cross-sections. In what follows, for the sake of brevity we 
present only results which include gluon shadowing. 

The nuclear ratio for incoherent scattering, — > jA', is shown in the Figure [3] It turns out that the effect of 
shadowing at small xb is twice as strong as for the coherent case. As a function of t, the nuclear cross-section is 
almost a constant-due to the nucleus breakup the nuclear cross-section is not suppressed by FAit) and decreases only 
as a function of the nucleon formfactor Fm (t) . The Q 2 dependence of the nuclear ratio is quite similar to what was 
found for the coherent case. 

For coherent RCS, jA — > jA, the results of evaluation are shown in Figure|4j Similar to DVCS, the cross-section is 
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Figure 2: The nucleus to nucleon cross section ratio for coherent DVCS as function of different kinematical variables. Left: 
shadowing as function of Bjorken xb with and without gluon shadowing for different nuclei, t = t m in, Q 2 = 5GeV 2 . Right: 
t-dependence for different A. x B = 10~ 3 , Q 2 = 5GeV 2 . 




Figure 3: Dependence of the nuclear ratio for incoherent DVCS on different kinematical variables. Left: IB-dependence, center: 
t-dependence, right: Q 2 -dependence. 



steeply falling with center of mass energy W (for DVCS xb oc 1/W 2 at fixed (Q 2 ,t)). As function of the momentum 
transfer t, the shadowing correction reveals the behavior qualitatively similar to the case of coherent DVCS (see Fig. [I} 
It steeply decreases similar to the behavior of the nuclear formfactor Fa (t) ■ 

For the incoherent RCS, 7 A — > jA', the results of evaluation are shown in Figured] Similar to the coherent case, 
the cross-section is decreasing as function of W, down to the values an order of magnitude smaller than would give a 
simple sum of the nucleon cross-sections. As function of t, the cross-section is almost constant because of absence of 
the nuclear formfactor. 



VIII. CONCLUSIONS AND PROSPECTS 



In this paper we considered DVCS and RCS on nuclear targets within the color dipole model. Results for the 
coherent shadowing are presented in Section IVIII We found that the magnitude of nuclear shadowing is large and is 
important for analysis of DVCS and RCS data. It may cause a substantial suppression of the nuclear DVCS cross- 
section. We observed that for incoherent case, the shadowing is stronger than for the coherent case. This happens 
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Figure 4: Nuclear ratio for coherent RCS as function of W (left) and t (right). 




Figure 5: Dependence of the incoherent nuclear RCS cross-section ratio on W (left) and t (right). 



since the incoherent cross section is proportional to the survival probability of the dipole, so it vanishes in the black 
disk limit (oc A 1 / 3 ). At the same time, the coherent cross section saturates at a value oc A 2 / 3 . For the RCS, the 
shadowing is larger than for the DVCS, since a real photon fluctuates to dipoles of a larger size than a virtual one. 

Our results for DVCS and RCS can be applied to eA collisions at future electron-ion colliders_(EIC and LHeC, [13, 
HH). Currently data for nuclear DVCS are available from the HERMES experiment at HERA [26|, however the values 
of Bjorken xb are too large to produce any sizable shadowing effects. Besides electron beams, the usual source of 
quasi-real photons, one can also use beams of charged hadrons. Provided that the transverse overlap of the colliding 
hadrons is small, i.e. the transverse distance b between the centers of the colliding particles exceeds the sum of their 
radii, b > R\ + R 2 , the electromagnetic interaction between colliding particles becomes the dominant mechanism. 
Such processes called ultra-peripheral collisions (UPC) can be studied in pp 7 pA and AA collisions. In particular, one 
can access RCS in the reaction 

At +A 2 ->• A 1 +j + A 2 . (42) 

The typical virtualities (Q 2 ,.) of the intermediate photon 7* are controlled by the formfactors of the colliding particles, 
and are small: 

where A is the atomic number of the hadron (A = 1 for a proton) interacting electromagnetically, i.e. emitting the 
photon, while the second hadron interacts strongly. 

Thus, (Q 2 ,*) is of the order of the soft hadronic scale, so the intermediate photon can be treated as a free Weizscker- 
Williams one, i.e. the amplitude of the process (|4"2")l can be described in terms of RCS. 
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These processes at the LHC will allow to study RCS at very high e nergies . The possibility of observation of such 
processes experimentally has been demonstrated by the STAR |60l - [62| and PHENIX [63| experiments at RHIC. It is 
expected that at LHC photon-proton collisions at energies up to y/s lp <8x 10 3 GeV can be observed (64|. 
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